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On the theory of superfluidity 
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We investigate the properties of dispersion spectra of one-dimensional periodic Bose systems with 
repulsive interparticle interactions. These systems with sufficient large interactions can support 
metastable supercurrent states, which correspond to the local minima of the dispersion spectra at 
non-zero momenta. The existence of local minima in the spectra and the energy barriers, which 
separate the minima, can be explained in terms of Bose exchange symmetry. We extend our study 
to the case of higher dimensional Bose systems. We suggest that superfluidity could be understood 
as a Bose exchange effect. 

PACS numbers: 03.65.-w, 05.30.Jp, 67.40.-w 



I. INTRODUCTION 

Superfluidity, a fascinating quantum phenomenon in 
which a system demonstrates frictionless motion, is one 
of the central interests in condensed matter physics. As 
a many-body quantum effect, however, an exact under- 
standing of superfluidity had remained elusive for more 
than six decades after its discovery in liquid ''He in 1938 
[l[ . It has be long suspected that superfluidity has some- 
thing to do with Bose exchange symmetry: the wavefunc- 
tion of a Bose system is invariant under the exchange of 
the coordinates of two particles (see, e.g., 0, 0, It 
has been pointed out that, especially in the case of one- 
dimensional Bose systems [3, fsj, if the many-body dis- 
persion spectrum of a system contains local minima at 
nonzero momenta, the many-body states corresponding 
these minima are metastable and current-carrying, i.e., 
they are supercurrents, since the energy barriers which 
separate the local minima prevent the supercurrents from 
decaying (see Fig. [T]for example). In this paper, we at- 
tempt to further argue that the existence of these minima 
of the dispersion spectrum is due to Bose exchange sym- 
metry. The paper is organized as following, in Sec. II 
we discuss a kinetic property of a one-dimensional Bose 
system with translational symmetry, and its connection 
with the critical velocity of superfluidity. In Sec. Ill we 
perform a perturbation analysis of dispersion spectrum of 
a Bose system in the weak interaction regime. In Sec. IV 
we compare the dispersion spectrum of a (spinless) Bose 
system with the dispersion spectrum of a Bose system in 
which the particles have some internal degrees of freedom 
(a spinor Bose system). With this comparison we illus- 
trate the role of exchange symmetry in determining the 
structure of dispersion spectrum of a strongly interacting 
Bose system. We briefly extend our study to the case of 
higher dimensional systems in Sec. V, while the detailed 
discussions on these systems will be presented elsewhere. 
The conclusion is given in Sec. VI. To avoid confusions, 
we restrictively define superfluidity as the property that 
a system supports persistent current in this paper. We 
don't attempt to link our theoretical model systems to 
any realistic systems, since our purpose is to seek general 



qualitative understanding of superfluidity. 



II. GALILEO INVARIANCE AND CRITICAL 
VELOCITY OF A SUPERFLUID 

We first consider N bosonic particles interacting with 
repulsive interaction, moving along a ring with radius of 
R, i.e., the system has a one-dimensional periodic geom- 
etry. The Hamiltonian has the form 
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where g > is the interaction strength, 9i is the angular 
coordinate, M is the mass of a particle and / describes 
the form of the interaction (/ > 0). We only consider 
short-ranged interaction and zero-ranged interaction. 

Three quantities, h, M and R, suffice to generate a set 
of units for the physical quantities involved in this paper, 
which will be used to render all quantities dimensionless. 
The unit of energy, however, is taken to be /2MR^ , 
where an extra factor of one half is involved for conve- 
nience. 

Due to the axial symmetry of the Hamiltonian, the to- 
tal angular momentum is a good quantum number. It 
is straightforward to show that the many-body spectra 
and eigen wavcfunctions, at two angular moments with 
one differing from the other by kN where k is an integer, 
can be related in the following way: giving an eigen wave- 
function -0(^1, 02, ■■■,&n) with angular momentum L' and 
with energy Etj, , a mapped wavefunction 
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will be an eigen wavefunction with angular momentum 

L' + kN (3) 

and with energy 

- E.^ + L^/N - L'^/N. (4) 



The factor in Eq. [2 f,ike,+ik92+...+ike« ^ describes 
center-of-mass motion and conserves the inner structure 
of ^/i in this one- dimensional periodic system. The energy 
difference between ^ and (p is independent of interaction, 
and one might attribute this kinetic property to Gahleo 
invariance (see also [1, 0] ) ■ 

We examine only the dispersion relation, E — E{L), 
where E is the lowest energy of the system at given L. We 
refer to E as the "yrast" energy, following the convention 
in nuclear physics. Generally, the yrast energy at i = is 
the global minimum of the dispersion spectrum, i.e., the 
ground state energy. Other local minima (if they exist) 
in the yrast energy at nonzero L, i.e., ... > E{L — 2) > 
E{L - 1) > E{L) < E{L + 1) < E{L + 2) < ... indicate 
metastable states carrying persistent currents [j, ^ . The 
peculiar possible existence of the local dispersion minima 
at none zero L is of quantum nature and beyond classical 
understanding. 

With the mapping (Eqs. [2]|4]), one will know the full 
dispersion spectrum if the part at angular regime < 
L' < N is giving (see Fig. [1]). For the energy difference 
between two neighboring yrast states, with Eq. [2] one 
easily obtains the following relations, 

E{L = L' + kN) - E{L - 1 = L' - 1 + kN) 

= E{L') - E{L' - 1) + 2k + (2L' -1)/N (5) 




FIG. 1: A schematic plot of many-body dispersion spectrum 
of a one-dimensional Bose system with repulsive interaction. 
Once the part of the spectrum at angular momentum regime 
< L < N (the thick line) is known, the rest of the spectrum 
can be obtained. The star marks the ground state, filled cir- 
cles (red or blue) mark the yrast states which can be mapped 
to the ground state, among which only those with angular 
momentum less than a certain value are supercurrent states 
(blue circles). See also 0]. 



and 

E{L ^kN+l)~ E{L -l = kN) 
= E{N + 1) - E{N) + 2k+l/N. (6) 

With large enough k, for example, k > max{E{j) — E{j + 
l))/2 (0 < j < N), the right hand sides of Eq. [5] and 
Eq. [6] are positive for all L. Hence, E{L — 1) < E{L) < 
E{L+1) < E{L+2) < ... < E{oo) with sufficient large L, 
the dispersion spectrum will not support any persistent 
currents further. This naturally explains that there is an 
upper limit of (angular) velocity of a superfluid, beyond 
which the system is not frictionless any more(see also 
Q). It will be shown that the yrast spectrum contains a 
local minimum at L = N, giving sufficiently large g, and 
has a local linear dispersion around this minimum. With 
the analysis above, it is clear that the slope of the linear 
spectrum at L = TV, or simply E{N — 1) — E(N), related 
to the slope of the spectrum at L = 0, will determine the 
critical (angular) velocity of the superfluid. 

Landau [7] ffist related superfluidity of a system with 
the properties of quasi-particle spectrum of the system. 
One shall distinguish Landau's quasi-particle spectrum 
and the many-body dispersion spectrum discussed above. 
A clear difference between two spectra can be seen if one 
compares their general behavior at rather large momen- 
tum p, where the quasiparticle spectrum likely approx- 
imates /2M but the dispersion spectrum is close to 

/2NM with N can be macroscopic (the ground state 
energy is shifted to be zero). 



III. THE STRUCTURE OF THE DISPERSION 
SPECTRUM AT WEAK INTERACTION REGIME 

For the purpose of simplicity and convenience, we shall 
also use a rotating frame of reference besides the labora- 
tory frame. In a rotating frame with angular velocity fJ, 
the Hamiltonian is — H — QL. 

A natural set of single particle orbits can be chosen to 
be ipm{0) = e™VV2^, m = 0, ±1, ±2, .... One easily sees 
that the yrast states of noninteracting system {g — 0) 
at the angular momentum regime < L' < N are the 
Fock states, \0'^1°),\0'^-^1^),\0'^-'^1^), ...,\0H^) (Here 
|0^°1^^) means that Nq particles occupy the orbit (po, 
and iVi particles the orbit (pi). Only two single-particle 
states ifQ and ipi are involved in these yrast states, while 
Fock states involving any other single particle orbits are 
excited states above the dispersion spectrum since they 
have relatively larger (kinetic) energy. In the rotating 
frame with Q — 1, ipo and pi are the degenerate single 
particle ground states with zero energy. 

In the limit of weak g, within the first order perturba- 
tion theory, the Fock states |0^-^'l^')(0 < L' < N) re- 
mains the yrast states. In the rotating frame with = 1, 
the kinetic energies of these Fock states are zero. The in- 
teraction energy of |0^~^ 1^ ), the same as the total en- 
ergy in the rotating frame, is calculated straightforwardly 
and can be written in four parts, 

E''^ = (N - L'){N - L' - l).ge/2 + L'{L' - l)ge/2 

+ iN-L')L'ge/2+{N-L')L'ge,j2 (7) 
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FIG. 2: Interaction diagrams for the Fock state |0^^^ 1^ ) . 
The dashed line represents the interaction, m is the angular 
momentum of the single particle orbit. In the last diagram, 
two particles exchange their orbits. 



with e = /^^ f{\9\)de, and with 
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One might think that the four terms on the right hand 
side (RHS) of Eq. [7]correspond to the four interaction di- 
agrams in the Fig. [21 respectively. For example, the first 
term counts the interaction energy with two (interacting) 
particles in the orbit with m = 0. Since there are N ~ L' 
particles in this orbit, out of which {N - L'){N - L' -l)/2 
pairs of particle can be formed. The sum of first three 
terms on the RHS of Eq. [7]is N{N - l)ge/2 and is inde- 
pendent of L' . The last term, causing by the exchange in- 
teraction, generates the parabolic energy barrier between 
the ground state and yrast state at L = TV (in the rotat- 
ing frame). Note that the argument here is equivalent 
to that of Nozieres in Q where he relates Bose-Einstein 
Condensation with exchange interaction. 

One might check with a system with = 2 to see how 
the exchange term originates. The wavefunction of the 
yrast state at L = is *l=o(6'i, ^'2) = e'°^ie*°^= /27r, and 
that at L = TV = 2 is *l=2(6'i, 6*2) = e^'^^e''^^ /2tt. Both 
wavefunctions have a product structure and are symmet- 
ric under exchange of two Bose particles (since two par- 
ticles occupy the same orbits). 

The wavefunction of the yrast states at L = 1, is 
*L=i(6'i,6'2) = \/2/2(e™ie^^= -h e**ie^o^^)/27r, which 
takes a form of a sum of two products, due to 
the symmetrization. The interaction energy, V — 
J J *L=i27i'5/(^i - d2)^L=idOid02 , second order in the 
wavefunction, will generate extra cross terms, given that 
the wavefunction itself is a sum of products instead of a 
product. 

If one ignores the symmetrization and considers an in- 
vahd wavefunction at L = 1, -^^^^f^d ^ e*0'^ie*''V27i", it 
will has the same interaction energy as that of ^'l=o,2- 
Comparing l=i a-nd the invalid wavefunction, one sees 
that symmetrization pushes up the (interaction) energy 
for the yrast state with L = 1. 



IV. A COMPARISON BETWEEN THE 
DISPERSION SPECTRUM OF A SPINLESS 
BOSE SYSTEM AND THAT OF A SPINOR BOSE 
SYSTEM 

In Sec. Ill, we show that the structure of dispersion 
spectrum of the Bose system is determined by the ex- 
change interaction in the limit of weak interaction. We 
shall show that the possible dispersion energy barrier be- 
tween the yrast state at L = and that at L = TV be- 
comes more apparent at larger and can be attributed to 
the exchange interaction. These are illustrated by com- 
paring the dispersion spectra of the Bose system consid- 
ered above, referred to a spinless Bose system, with the 
dispersion spectra of a spinor Bose system, which only 
differs from the spinless system by allowing two degener- 
ate internal states |a) and |/3) to each particle. 

We can formally introduce (pseudo-) spin operators 
, Sz associated with the internal states. For example, 
Sz\a) = |a)/2 and s^,]/?) = -|/?)/2. We then refer to \a) 
as spin up and \(3) as spin down. The interaction does 
not affect the spin and the spinor Bose system has SU (2) 
symmetry. In second quantization, the Hamiltonian is, 

l,(y=a.,P 

k^l,rn^n,(T,ij' 

where cj ^ ( ci^a) is the operator which creates (annihi- 
lates) a particle at the single particle state (pi{9)\a), and 
where 
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Due to the SU{2) symmetry, S'^ = (si + S2 + ... + sn)^ 
and Sz = Szi + Sz2 + --- + SzN are good quantum numbers. 

Let us first examine how the spin degrees of freedom 
interplay with the exchange interaction. A many-body 
wavefunction of a Bose system is required to be symmet- 
ric under the exchange of any two particles. If the wave- 
function is symmetric/antisymmetric under the exchange 
of their spin degrees of freedom, it also has to be sym- 
metric/antisymmetric under the exchange of their spatial 
degrees of freedom. For example, we consider a system of 
two Bose particles which occupy different single-particle 
spatial orbitals, ipk{9) and Lpk+i{9). The wavefunction 
with S* = is 

1^ = 0,5, =0) = ibfc(0l)(^fc+l(02)-95fc(02)(/?fc+l(0l)] 

x(|ai)|/32)-|ai)|/32)I^l) 
while the wavefunction with 5 = 1 and = is 



\s = 1,5, = 0) = -[ipk{ei)^k+i{e2) + Me2)^k+i{ei)] 



x(|ai)|/32) + K)|/?2)II2) 
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The state with 5* = has an interaction energy of 
ge — ge^x since the exchange energy is equal to —ge^x- 
For the state with 5 = 1, the exchange energy is gcex and 
the total interaction energy is ge + geex- In this system, 
the exchange interaction favors the S* = state, which 
can be understood easily. With 5 = 0, the spatial part of 
the wavefunction is antisymmetric and it vanishes with 
d\ = 02, implying that one particle has a small chance 
to approach and to interact effectively with the other 
particle. 

To demonstrate the spin structure of the TV-particle 
yrast states of the spinor Bose system, we will once again 
discuss the perturbative regime first. At small g and a 
given L', in first order perturbation the spatial orbital 
configuration consists oi N — L' particles occupying the 
ipo state and L' particles occupying the ipi state. Accord- 
ing to the simple argument given above, it is energetically 
favorable for the yrast state to accommodate as many 
spatially antisymmetric pairs as possible. For L' < N/2, 
the maximum number of such pairs is L' and the rest 
N — 2L' particles are in the state (/Sq- Each pair has 
a zero contribution to the total spin and the unpaired 
particles are symmetric in spin, thus the total spin is 
{N-2L')/2. For N/2 < L' < N, S = {2L' - N)/2. Our 
numerical calculation confirmed this spin structure. In 
addition, we found for the yrast spectrum as a function 
of L', = —(TV — 2S)geex, apart from a constant term 
N{N — l)ge/2 (This constant term will be ignored there- 
after). Here takes the form of an exchange energy, 
where each spatially antisymmetric pair contribute with 
an exchange energy of —gtex ■ 

We shall consider the spinor system with any g. It is 
found numerically that the total spin 5 of a yrast state 
at a given angular momentum L' , as a function of g, is 
stable, i.e., it does not change for any finite g. It might 
be understood in the following way. The exchange in- 
teraction plays a dominant role in determining the low 
energy spectrum. At large g, even though the yrast state 
involves more single-particle spatial orbits and all orbits 
contribute to the exchange energy, the leading contribu- 
tion still comes from tpo and ipi. To reduce the leading 
contribution, S will therefore remain the same as in the 
perturbative regime. 

For the spinor Bose system, E^{L'), is linear in g for 
small g, then decreases monotonically with g and finally 
saturates to —L'{N — L')/N at sufficient large g, as we 
found numerically [9|. Correspondingly, the dispersion 
relation does not have any local metastable minima in 
the laboratory frame. Hence, the spinor Bose system 
does not support persistent currents [10]. 

We shall compare the spinless Bose system and the 
spinor Bose system. One can see that a many-body state 
of the spinor system with S = Sz = N/2 is equivalent 
to a many-body state of the spinless system, since all 
spins point to the same direction and the spinor system is 
equivalently spinless. The yrast spectrum for the spinless 
system corresponds to the lowest energy eigenstate of 
the spinor Bose system at a given angular momentum. 
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FIG. 3: The yrast spectra in the lab frame. We use dots 
for the spinless Bose system and crosses for the spinor Bose 
system. Left panel refers to a contact interaction and right 
panel to a finite range interaction .f(\0\) = e~* ^'^'^ /V 47re with 
e = 0.05, illustrated in the inset figure. In both graphs, A'^ = 
10 and g — 1. 
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FIG. 4: The yrast spectra in the rotating frame, relative 
to the yrast energy at L = 0. Left panel refers to the 
delta interaction and right panel to a finite range interac- 
tion f{\e\) = e-''^/'^'/V^ with e = 0.05. The total spin 5" 
of yrast states is written out in the left panel (5 = N/2 for 
spinless system). In both graphs, N — 10 and g — 1. 



with S = N/2 and Sz = N/2. The exchange interaction 
associated with S = N/2 pushes up the energy of yrast 
state of the spinless system, relatively to the energy of 
the yrast state of the spinor system which has S = \N — 
2L'\/2. This generates a barrier between the yrast state 
at L' = and that aX L' — N in the rotating frame. 
For g larger than a critical value, the yrast spectrum of 
the spinless Bose system develops a minimum sA L = N 
in the laboratory frame, allowing for the possibility of 
metastable persistent currents. 

As a manifestation of the exchange interaction, for a 
spinless Bose system E^ increases with the (exchange) 
interaction strength g. This motivates the following re- 
lation, 



E'\L',g2)>E'\L',gi) if g^ > gi 



(13) 



At large iV and small angular momentum L' = Z > 0, we 
assume that E^ has a power series expansion as follows. 



E^{l,g) = ao + ai{g)l/N + a2{g){l/Nf 



(14) 
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FIG. 5: Critical angular velocity uoc. as a function of g. The 
dashed line is for the delta interaction, and the solid line for a 
finite range interaction f{\0\) = ''"''^/\/47re with e = 0.05. 
In both cases N = 10. 



ao = since E^{1 = Q,g) = Q (the ground state energy 
is shifted to be zero). In the perturbative regime (small 
g) with E^{l,g) = 1{N - l)gt^^ we find ai{g) = gN'^eex 
which is not zero. Therefore the linear term is the lead- 
ing (nonvanishing) term of E^{l,g). ai{g)/N might 
be identified with critical angular velocity lUc, and thus 
LUc = gNe^x in the perturbative regime, we ignore all 
higher terms other than the linear term at the RHS of 
Eq. [M] at small l/N. Combining the two relations [T3l 
[Til we find 

ai{g2)/N > ai{g^)/N if 52 > 5i (15) 



Wc(ff2) > Wc(ffi) if 92 > 91- (16) 

Therefore, the critical angular velocity, as a function of 
g, is linear at small g and increases monotonically with 
g. In order for the system to support a supercurrent at 
L = N , LUc needs to be large than unity. 

Some numerical examples are plotted in Fig. [3] and 
Fig. m The yrast spectrum of a spinless system and 
that of the corresponding spinor system are presented in 
left panel of Fig. [31 For both systems, N = 10, f{\9i - 
^2!) — S{9i — O2) (,i.e., a contact interaction), and g = 
1 {Ngeex = 10 ^ 1, and therefore the interaction is 
beyond the perturbative regime). It is clear that there 
is a supercurrent at L = 10 for a spinless system while 
there is no supercurrents for the corresponding spinor 
Bose system. For both systems the corresponding E^^ is 
presented in the left panel in Fig. [4l The right panels of 
Fig. [21 and Fig. [H correspond to the systems with finite 
range interactions, f{\0\) = ^'^'^/y/Aire with e = 0.05. 
For both systems, = 10, g = 1. 

We also calculated the critical velocity ojc of spinless 
systems and plotted Wc vs g in Fig. \E\ 
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FIG. 6: The geometry of a 2D periodic system. The edge 
AB is identified with DC and BC identified with AD. The 
linear size of the system is 2-kR. 



V. HIGHER DIMENSIONAL SYSTEMS 

The above understandings of dispersion spectra of one- 
dimensional Bose systems can be generalized to the case 
of higher dimensional systems straightforwardly. Here 
we just briefly discuss 2D periodic systems (see Fig. [6]) 
with finite range interactions while the more detailed dis- 
cussions and the numerical checks on higher dimensional 
systems will be presented elsewhere. 

The total momentum components P^, Py of a 2D sys- 
tem are good quantum numbers and we shall consider 
only the dispersion spectrum E — E{Px, Py) where E is 
the lowest eigen level at a given momentum. Local min- 
ima (if they exist) of the dispersion spectrum at non-zero 
momenta will indicate the supercurrent states. 

The perturbative analysis of the dispersion spec- 
tra in the limit of weak interaction strength g is 
straightforward, similar to the case of one-dimensional 
systems. For example, the yrast states of a spin- 
less system at the momenta < Px < N,Py = Q 
will be approximated well by Fock states 
1(0,0)^(1,0)0), 1(0,0)^-1(1,0)1), 1(0,0)^-2(1,0)2), 

1(0, 0)°(1, 0)^), where {m,n), a pair of inte- 
gers, denotes the single particle orbit ipm,n{x,y) = 
l/(27r)e*"^/2'^'^+™i'/27r_R^ rpj^g interaction energies of 
these states, as a function of Px will has a parabolic part 
(the exchange term) c>c Px{N — Px) besides a constant 
part (the Hatree term). If g is above a critical value, 
the dispersion spectrum will contain a local minimum at 
Px — N,Py = 0. By contrast, the dispersion spectrum 
of a spinor Bose system will contain no local minimum 
other than the global one at Px = 0,Py = 0. 

Due to the Galileo invariance, the spectrum at large 
momenta can be written out by mapping them to the 
spectrum at momenta regime < P^; < and < 
Py < N. Following the mapping algebra ( 2D versions 
of Eqs. 7-9), one can find that the dispersion spectrum 
of a spinless Bose system will contain no local minimum 
at given large enough momentum, i.e., there is an upper 
velocity limit for a superfluid. 

One rather new finding about the 2D systems is fol- 
lowing: for a spinless system with A'^ = where j is 
an integer and with given large enough g, there will be 
also supercurrent states at Px = kxj,Py = kyj, where 
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kx^ky = 1,2, 3,..., J — 1. One might note that, with 
N = j'^ , the system allows to form a lattice (crystal- 
lization) at large repulsive interactions. The new super- 
current states, written as superpositions of many-body 
Fock states, can not involve any BEC-like states (Fock 
states with all N particles occupy the same orbit). The 
origin of the new supercurrent states, can be understood 
in terms of Bose exchange symmetry (the explanation is 
presented elsewhere). 

VI. CONCLUSION 

We have studied the properties of dispersion spectra 
of Bose systems with repulsive interactions. We argued 
that the possible existence of local minima of the spec- 
tra and the energy barriers separating these minima are 



due to the Bose exchange symmetry. Although the nu- 
merical checks can be only done to systems with limited 
N, the microscopic physical arguments can be applied to 
systems with any large N. 

The quantum phenomenon of superfluidity is rather 
puzzling. However, one might get a simple picture of su- 
perfluidity if one, i) relates it with dispersion spectra of 
the superfluids 0, Q, ii) explains the non-classical fea- 
tures of the dispersion spectra in terms of the exchange 
properties of Bose particles. One might draw an analog 
between this picture of superfluidity and Heisenberg's ex- 
change theory of ferromagnetism In fact, part of 
arguments presented in Sec. IV bears a large similarity, 
both formally and in content, to Heisenberg's arguments 
on ferromagnetism. In conclusion, we suggest that super- 
fluidity might be understood as a Bose exchange effect. 
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